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Abstract 

We study the cosmological solutions of the one loop corrected superstring ef- 
fective action, in a Friedmann-Robertson -Walker background, and in the presence 
of the dilaton and modulus fields. A particularly interesting class of solutions is 
found which avoid the initial singularity and are consistent with the perturbative 
treatment of the efi'ective action. 



* Research supported by EEC contract Ref B/SC1*-915053. 
** Laboratoire Propre du Centre National de la Recherche Scientifique UPR A. 0014. 



1. Introduction 



Einstein's theory has been very successful as a classical theory of gravitational 
interactions. However, a quantum theory of gravity would require to incorporate 
Einstein's theory in a more general framework. Today our best candidate for such 
a framework is superstring theory which also has the prospect of unification of 
all other interactions. Superstrings appear to involve a minimum length which is 
of the order of the Planck scale and, thus, they are expected to lead to drastic 
modifications of the Einstein action at short distances. These modifications can in 
principle have cosmological consequences which distinguish string cosmology from 
the standard model and provide some indications for a stringy origin of the universe. 
Furthermore, one hopes to find resolutions to some puzzles in Einstein cosmology 
including the initial singularity problem. 

String theory gives rise to two kinds of such modifications. The first is associ- 
ated with the contribution of the infinite tower of massive string modes and leads to 
Qj'-corrections, while the second is due to quantum loop effects. Both these contri- 
butions can be studied in the context of an effective field theory which involves only 
the massless string modes. The effective Langrangian can be derived from string 
theory using a perturbative approach in both the string tension a! and string cou- 
pling expansion. An alternative direct way to take into account all a'-corrections 
at the classical string level is to consider conformal field theories which describe 
exact string solutions in time-dependent gravitational backgrounds [1,2]. However 
this is not in general sufficient for probing the short distance behavior of the theory 
since the gravitational coupling is dimensionful and quantum corrections become 
important in this region. 

One of the unique properties of the string effective action is that couplings are 
field-dcpcndcnt. For our purposes, the relevant fields are the dilaton which plays 
the role of the string loop expansion parameter, and the moduli whose vacuum 
expectation values describe the size and the shape of the internal compactification 
manifold. Since these fields have no potential their contribution is expected to be 
important in any cosmological solution of the effective action. Here we restrict 
ourselves to the simple case of a single modulus field, besides the dilaton, which 
corresponds to the common compactification radius. 

The tree-level string effective action has been calculated up to several orders in 
the a'-expansion in both the sigma model approach [3] , where one considers strings 
propagating in background fields, and the S-matrix approach, where the effective 
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action is computed directly from string scattering amplitudes [4] . It turns out that 
there is no moduli dependence of the tree-lcvcl couplings. The one loop corrections 
to the gravitational couplings have been also calculated recently in the context of 
orbifold compactifications of the heterotic superstring [5]. It has been shown that 
there are no moduli dependent corrections to the Einstein term while there are non- 
trivial contributions. They appear as the Gauss-Bonnet combination multiplied 
by a function of the modulus field. It is interesting to consider the implications 
of this term to the the cosmological solutions of Einstein equations for two main 
reasons. Firstly this term is subject to a non-renormalization theorem which implies 
that all higher loop moduli dependent contributions vanish. On the other hand, 
it breaks the continuous isometries of the tree level modulus kinetic terms leaving 
intact only the duality symmetries. 

The purpose of this paper is to study the evolution of the equations of motion of 
the corresponding effective Lagrangian. We investigate the asymptotic solutions and 
we show that the one loop moduli dependent action contains a class of interesting 
cosmological solutions which avoid the initial singularity. This is possible only 
for a definite sign of the corresponding four-dimensional trace anomaly for which 
the strong energy conditions [6] related to the modulus energy-momentum tensor 
can be violated. These solutions start from flat space-time in the infinite past, 
they pass through an inflationary period and they end up as a slowly expanding 
universe. Although in our analysis we omit terms higher than quadratic in the 
Riemann tensor, we will argue that the above result persists in the full theory 
under certain assumptions on the moduli dependence of loop corrections to higher 
derivative terms. 

The Gauss-Bonnet density multiplied by the dilaton field is already present at 
the string trcc-lcvcl in the next-to-lcading order of the a'-expansion and it has been 
studied extensively in the literature [7]. As we will show in Section 3, the addition 
of this term alone does not lead to violation of the energy conditions and thus it 
cannot provide any singularity-free solution, in contrast to the modulus dependent 
loop correction. 

We also examine the possible existence of time dependent solutions which fix 
asymptotically the vacuum expectation value of the modulus field. This is moti- 
vated from the fact that the Gauss-Bonnet contributions could be viewed as a time 
dependent potential for the modulus with an extremum at the self-dual point. We 
show that this possibility can be realized only at early times and in singular strong 
coupling solutions. 
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This paper is organized as follows. In Section 2 we review the string effective 
action and the loop corrections to the gravitational couplings and wc derive the 
equations of motion for the coupled system of graviton, dilaton and modulus field. In 
Section 3 we study the cosmological solutions in the simple case where the dilaton is 
ignored. We classify all asymptotic solutions and we show by numerical integration 
that two of them can be smoothly joined avoiding the singularity. We justify this 
behavior by demonstrating the violation of the weak and strong energy conditions 
related to the singularity theorems. Finally, we explore the parameter space of all 
initial conditions by the use of the corresponding phase diagram. In Section 4 we 
repeat the same analysis in the presence of the dilaton field and show that the 
previous results remain unafi^ected. Our conclusions are summarized in Section 5. 



2. The loop corrected effective action and the equations of motion 

Let us consider the universal part of the effective action of any four-dimensional 
heterotic superstring model which describes the dynamics of graviton, dilaton S and 
the common modulus field T. At the string tree-level, and up to first order in the 
a'-expansion, it takes the form [3,4]: 

1 D'lD^ DTDT 1 1 ~ 

where D denotes covariant differentiation, R is the scalar curvature, k = ^/St^Gn 
with Gjsf the Newton's constant, Rq^ is the Gauss-Bonnet integrand 

and RR = V^^^^R,.^''^ RnXar-^ The inverse of Res' plays the role of the string 
coupling constant squared while ImS is a pseudoscalar axion. Finally the real part 
of the complex modulus field T corresponds to the square of the compactification 
radius. 

At the one loop level the moduli dependence of the gravitational couplings 
in the case of the heterotic string compactified on a symmetric orbifold has been 



1 We use the conventions i?^^, = i?^^^ , R^^^ = dj"^^ + . . ., and r/"^^"^ 
with e^^^'' = -eijk. 
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studied in [5]. It is shown that there are no moduh dependent corrections to the 
Einstein term, while the contributions to the four derivative gravitational terms 
take the form 

A£eff = A(T, f)R%B + Q{T, f)RR . (2.3) 
The moduli dependent functions are defined as 

^iT,T) = ^ln[{T + f)\rj{iT)n (2-4) 

and G(T,f ) = -zA(T,f ) , where //(r) = q^/^'^Un>i (l " ?^'')' with q = e*'^^, is 
the Dedekind ry-function. A(T, T) is invariant under the duality SL{2, 7i) transfor- 
mations T — > 1/T and T ^ T + i which are the discrete subgroup of the continuous 
5'L(2,R) isometry group of the modulus kinetic terms in (2.1). The coefficient bgr 
is proportional to the four dimensional trace anomaly of the N=2 sectors of the 
theory 

i(-3iVv + Ns) - y (-3 + , (2.5) 

where Ns, Ny and A^3/2 denote the number of chiral, vector and spin-3/2 massless 
supermultiplets . 

We now consider the effective action (2.1)+(2.3) in a spatially flat homogeneous 
and isotropic Robertson- Walker background. Since RR vanishes identically in this 
background, it is consistent with the equations of motion to assume for simplicity 
ImS = constant and ImT = 0. The latter is required in order to have vanishing 
derivative of A(T, T) with respect to ImT. Setting the length scale k = 1 and 
defining KeS = p-e'^, where g is the four-dimensional string coupling, and ReT = 
e^'^ the effective action takes the form 

R l,„^,o 3,„,o 1 



5 = d^x 



^ + :,{D^r + -ADaf + -(A - 6 i{a))Rt,^ 



(2.6) 



where A = ^, 5 = |^ and ^(a) = In [2e''r]'^{ie 

The equations of motion derived from the action (2.6) by varying with respect 
to the metric, the dilaton and the modulus field have the form 

^ -lD'D,a+^4-Rl^ = ^, (2.8) 



-g5(7 2 da 
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where 



(2.9) 



(2.10) 

and/ = 3^(A e'^-5i[a)). 

Substituting the spatiaUy flat Robertson- Walker ansatz for the metric 

g^,^{l,-^^5i,) (2.11) 

in (2.7)-(2.10) and considering only time dependent fields, we end up with the 
equations 

0-2 $2 

3a;^ - 3— - — + 2Afu^ = , (2.12) 

6^ • •• 

2cl; + 3a;2 + 3— + — + 16/0;^ + 8fco^ + IGfuiu = , (2.13) 

a + 3uj& + 5^ib^{uj^ + u) = , (2.14) 
oa 

^ + 3cj6-3Ae*cj^(a;2 + a)) = , (2.15) 

where (2.12) and (2.13) correspond to the Tqo and components of (2.7), re- 
spectively. These two equations are not functionally independent because of the 
Bianchi identity related to the conservation of the total energy momentum tensor. 
In fact the linear combination a;x(2.13)-d;x(2.12) -|(7x(2.14)-^6x(2.15) yields 
the time derivative of (2.12). Thus, wc can reject (2.13) and consider the system of 
(2.12), (2.14) and (2.15) as the independent equations of motion. 

Finally, let us also define here the energy density (p) and pressure (p) of the 
dilaton-modulus matter system. Assuming a perfect fluid form for their energy 
momentum tensor Too = p, Ta = pe^^ ^ and using (2.12)-(2.13), we get 

p =3a;^ 

(2.16) 

p = - (2a; + Zu^) . 
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3. Analysis of the metric-modulus system 



For simplicity we start our analysis by neglecting all the dilaton related terms. 
Although $ = constant is not a solution of the dilaton equation of motion (2.15), the 
metric-modulus system provides a simple model to study the effect of the Gauss- 
Bonnet loop correction (2.3). This simplification will be justified in Section 4, 
where the full system will be examined, and we will show that the main results of 
the following analysis remain valid even in the presence of the dilaton. 

The system of equations of motion (2.12)-(2.15) is now reduced to 

Au^ -&^- 25^&u;^ = , (3.1) 
da 

a + 3u;& + S^dj^(u;^ + C0) = . (3.2) 
oa 

Note that the absolute value of the trace anomaly coefficient 5 can be absorbed by 
a time rescaling 

t^t' = ^/\5\ t (3.3) 

which implies that we can replace 6 by its sign. From the expression (2.5) one sees 
that in any theory with N < 4 supersymmetry the sign of S tends to be positive 
unless there is a considerable excess of vector bosons. 

The non-linear system (3.1)-(3.2) can be integrated using numerical methods. 
Given initial values for cr, uj and uj, the value of & can be determined from (3.1). 
Then, the time derivative of (3.1) together with (3.2) form a second order system 
of differential equations which can be numerically solved given the above boundary 
conditions. 

Before proceeding to the numerical integration, we shall first derive analytically 
the asymptotic solutions in the limits t ^ oo and t ^ 0. For this purpose we will 
use the asymptotic expansion of ^'{(j) = |^ for cr — > ±oo: 

e'(a)~-sign(a)|el-l . (3.4) 

In the limit t — > oo we find two kind of solutions. The first is obtained using the 
ansatz 

CO = ujn + a\nt , 

(3.5) 

(J = ao + pint , 
with uq and ctq constants, and it leads to two possibilities: 

(Aoo) : a=^,m = ^ (3.6) 
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and 



(Boo) : \/3\=2, a^-a^ + 5a-l = 0, e''^^'^^''' = ^^j^^ , (3-7) 

which leads to 5 > , a ~ 0.207 and sign(/?)a"o ~ 4.64 — In 5. (Aqo) is actually 
the asymptotic solution of the tree level system (5 = 0), while (Boo) is a new 
asymptotic solution where the Gauss-Bonnet term is important at late times. They 
both describe a slowly expanding universe, while the radius e^"^ is either expanding 
or contracting due to the duality symmetry a — > —a of equations (3.1), (3.2). 
The second kind of solutions is obtained using the flat space ansatz 

U> = Uq + UJit"' , 

(3.8) 

a = ao + pint , 

which leads to 

1 '^ 

(Coo) : a = -l,m = 5, e^'^^^f'^o = (3.9) 

This describes an asymptotically flat universe with slowly expanding (or con- 
tracting) radius. 

In the limit t — > we flnd only one asymptotic solution: 

CO —loq -|- a In t 

(Ao) : (3.10) 

cr =(To + ait^ 

with 

a=l, P = 2, = — (3.11) 

This is a singular solution which also fixes the modulus field to a constant (but 
arbitrary) value. 

Note that the modulus equation of motion (3.2) for slowly varying cu can be 
considered as describing the motion of a in the presence of a potential proportional 
to ^(cr). This potential has an extremum at the self-dual point a = which corre- 
sponds to a minimum for 5 < 0, providing a possible mechanism to fix the value of 
the compactification radius. Unfortunately, the first equation (3.1) requires u> — 
when cr = which leads to vanishing potential for a. The resulting solution a = 
and (jj = constant cannot be continuously approached in the asymptotic region. In 
the next Section we will see that in the presence of the dilaton, (3.1) has a slowly 
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expanding solution {uj Int) which leads to a realization of this mechanism at 
early times. In addition, the possibility of fixing the modulus through the same 
mechanism at finite time, in a region where cb exhibits an extremum, remains open. 

Numerical integration of the system (3.1)-(3.2) verifies the existence of the 
above list of asymptotic solutions but it also reveals another very interesting char- 
acteristic : For 5 < 0, there exists a region of boundary conditions, for which the 
two asymptotic solutions A^o and Coo are smoothly joined avoiding the singularity. 
In fact, starting from the asymptotically flat solution Coo at the inflnite past one 
is always continuously driven to the slowly expanding solution Aoo at the inflnite 
future. For the rest of boundary conditions, as well as for 6 > 0, the singular so- 
lution (Aq) is recovered. A typical non-singular solution is presented in flgs. 1-2? 
Fig. 1 shows that the expansion rate of the universe u starts from a zero value (flat 
space time) at t — > — oo, it grows up to a maximum value, and then it falls down 
again as 1/t at t — > oo. The scale factor, correspondingly, starts from a constant 
value, it goes through a period of rapid expansion (inflation) and it ends up to a 
slowly expanding universe. On the other hand, flg. 2 shows that the modulus fleld 
starts from — oo corresponding to zero compactiflcation radius at the remote past, 
it passes through the self-dual point a = during the inflationary period, and it 
ends up to a slowly expanding regime at the inflnite future^. The vanishing of a in 
the region where iv is maximum is a general feature of this kind of solution, which 
is due to the form of the Gauss-Bonnet loop correction as described above. 

Figs. 1,2 also show that the obtained non-singular solutions can have all time 
derivatives of uj and a less than one in Planck units which is consistent with our 
approximation of neglecting higher derivative terms in the effective action. This is 
in general sufficient provided that the moduli-dependent coefficients of the higher 
order terms are not large enough to compensate the derivative suppression. Under 
this assumption the main features of these solutions are expected to survive after 
higher order corrections are taken into account. 

The avoidance of the singularity is accompanied, as expected, by a violation of 
the weak and strong energy conditions [6] related to the energy-momentum tensor of 
the modulus fleld, which are illustrated in flg. 3. This violation can be demonstrated 
analytically in the following way. The system of equations consisting of the time 
derivative of (3.1) together with (3.2) can be solved for u and then ^'(cr) terms can 



^ In all plots we have used \5\ = 64. 

^ The duality symmetry implies also the existence of the dual solution a — > —a. 
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be removed using (3.1). Using (2.16) we obtain 

16a;* — Sw^cr^ + ocr* -1 ^\ 

p + 3p = -6(.; + ^) = ^ 16c0^_ 8^2^2 + 5^4 

where ^"{cr) = g2|:- It is easy to see that since ^"{cr) < 0, equation (3.12) imphes 
that p + p > and p + 3p > for 5 > 0. Thus, the energy conditions can never 
be violated in this case and we cannot avoid the initial singularity. This is also the 
result of the numerical integration, for S > 0, which always leads to the singular 
solution Ao at early times. 

On the contrary, as seen from (3.12), 6 < allows for solutions which explicitly 
violate both energy conditions. Of course this is in general a necessary but not 
sufficient requirement in order to avoid singularities. In our case, the numerical 
integration of the system (3.1)-(3.2) has shown that the violation of the energy 
conditions at some instant is also sufficient. In fact, the non-singular solutions 
can be obtained by imposing a) ~ p + p = 0, which implies p + 3p < 0, at the 
starting point of our integration. The initial values for a and cb at this point obey 
a constraint which can be derived using equations (3.12) and (3.1). For a specific 
value of the modulus a the expansion rate uj is given by 

4zy^2^ - biy^z^ + 40iy^z^ - 8(2 + 7iy'^)z'^ + 96^ - 144 = , z>^ (3.13) 

where z = (^C"(^)'^max u = ^rr^ satisfying —iKuKl due to the properties of 
^(cr). One can show that the above equation has at least one solution for z for any 
value of v. In terms of our original variables a and u>, these solutions are plotted in 
fig. 4. 

The entire region of initial conditions which lead to singular or non-singular 
solutions for 5 < can be explored by the use of the corresponding phase diagram. 
In fact, the system of second order differential equations (3.1)-(3.2) can be reduced 
to a single first order equation for the variable cj as a function of a. First (3.1) can 
be solved for & in terms of uj and a: 



a = -6C'{a)u^ ± \u;\y/6^$'{a)^u;^ + 4: . (3.14) 

On the other hand, the time-derivative of (3.1) together with (3.2) leads to (3.12) 
which can written in the form: 



da & IGu^ - 8u^&'^ + 5d4 



(3.15) 
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Inspection of (3.14) shows that the two branches defined by the two signs of the 
square root correspond to two disconnected classes of solutions associated with 
positive or negative values of a. These are related by a duality transformation 
o —>■ —a, and thus we can restrict to the case d > 0. 

The phase diagram of (3.15) is presented in fig. 5. The asymptotic solutions 
Coo and Aoo are located in the regions a — oo and a +oo, respectively, while 
the singular Aq appear in the region ci; — oo, cr > 0. As one sees, for negative 
initial values of a non-singular solutions are obtained for all values of lo. They 
always interpolate between an asymptotically flat and a slowly expanding universe. 
The maxima of oj correspond to points of the curve {& > 0) shown in fig. 4. On 
the other hand, for positive initial values of a singular solutions are obtained for all 
values of uj lying above a critical curve (bold line). Note that they never end up to 
asymptotically flat space time. Comparing to other known non-singular solutions 
[8,9] derived from higher derivative effective gravity theories, our solutions extend 
over an infinite region of the phase space where the maximum value of curvature is 
not bounded. 

As already mentioned, another interesting characteristic of these non-singular 
solutions is that they contain an infiationary regime.^ This is expected at least 
during the period when the energy conditions arc violated due to the development 
of negative pressure. The numerical analysis shows that the amount of infiation 
increases with rising values of ^^/^c^max given in fig. 4. For values of order one 
only a few e- foldings are obtained (see fig. 1), while larger values of this parameter 
break the validity of the perturbative treatment of the effective action. In any case, 
these solutions provide an example of a cosmological model which is driven to an 
inflationary era and then exits in a flnite time. This is in contrast to the behavior 
of other non-singular solutions which start from de Sitter space in the remote past 
[8,9]. In fact the solutions presented here start asymptotically from flat space-time 
which is welcome as it has been recently argued that eternal de-Sitter inflation may 
not be possible without a beginning [11]. Moreover, their characteristics are very 
similar to those arising in the context of a "pre-big-bang" scenario motivated by 
generalized scale- factor duality symmetries of string theory [2]. 

The properties of the non-singular solutions depend crucially on the form of 
the modulus dependent function ^(cr) arising from the string loop corrections to 

^ A discussion of the contribution of the string loop moduli dependent corrections 
to the gravitational couplings in connection with the inflationary solutions of ref. 
[8] was reported in ref. [10]. 
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the dimensionless gravitational couplings (2.6). For instance if ^(cr) ~ o"^ de Sitter 
space is obtained asymptotically. On the other hand if one considers the dilaton 
field instead of the modulus in the presence of the first order in a' Gauss-Bonnet 
interaction (2.6), the singularity can never be avoided. In fact the analysis is equiva- 
lent to that of the metric-modulus system (3.1)-(3.2) with the substitutions cr — ^ 
^C'(^) ~^ — Ae*^. Then, the energy conditions (3.12) have also the same form with 
— 5^"(cr) replaced by Ae* and they cannot be violated since \ — is positive. 



4. Analysis including the dilaton terms 

We now extend our analysis including the dilaton contributions in the effective 
action (2.6), which leads to the equations of motion (2.12), (2.14) and (2.15). A 
time dilatation combined with a shift in the dilaton field, 

$ ^ = ^ + 2\n\S\ - In A 

can be used to eliminate both A and the absolute value of 6. In analogy with 
the analysis of the previous Section we shall first derive the asymptotic solutions. 
As one can see from (3.4), the asymptotic behavior of the modulus dependent 
Gauss-Bonnet coefficient ^{a) for large a is proportional to el'^l, while the dilaton 
dependent coefficient is proportional to e* . It follows that the asymptotic solutions 
of the metric-modulus system derived in Section 3 will survive in this case, with the 
dilaton either being negligible or behaving similarly to the modulus. In addition, 
some new solutions are expected for dominant large asymptotic values of the dilaton 
field. 

In the limit t ^ oo the asymptotic solutions (Aqo) and (Boo) of Section 3 are 
extended with the dilaton being of the form $ = 0o + 7lnt to the solutions: 

{A'^):a=^, 9/32 + 37^=4 (4.2) 

and 

(B'^) : = 7 = 2 , a ~ 0.205 , sign(/3)cTo ~ 4.67 -ln6 , 6>0 , 0o ~ 3.62 . 

(4.3) 
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Furthermore, (Aqq) can also be extended with the dilaton of the form $ = (f)o + (f)it'^: 

(A'4): a=i, |/?| = ^, 7=-2, 01 = -^. (4.4) 

In the above solutions the dilaton field either grows to plus or minus infinity (A^ 
or and B'^) corresponding to weak or strong string coupling, respectively, or it 
reaches asymptotically an arbitrary constant value (A'^). As in the metric-modulus 
case in the solution A'^ the Gauss-Bonnet terms become irrelevant at large times. 

Similarly, the asymptotically flat solution Cqo can be extended with the dilaton 
being either of the form $ = (pQ + (ffif^: 

{C'^) : a = -1 , |/?| = -7 = 5 , = ^e^^o;? , e-sn(/3> 

o 

or of the form ^ = (pQ + 'ylnt: 
{C':,) : a = -1 , = 7 = 5 , 5 > , e-^"(/')-° = , e^o = _1_ . (4.5) 

In the dilaton goes to a constant, while C'^ is a weak coupling solution. 

There are also three new solutions where the modulus field goes asymptotically 
to an arbitrary constant. They are obtained from A'^ (Dqo, D'oo) ^'00 (Eoo) by 
interchanging the role of a and ^, a — ao + ait^ and ^ — (po + 'jlnt : 

(Doo) : a = ^ , |7| = ^ , /3 = -2 , ai = j^^^'l^o) , (4.7) 

(D;„) : a ~ 0.223 , 7 = 2 , /3 = -2 , (Ji - 0.0025^'(cto) , (po ~ 3.24 , (4.8) 

(Eoo): a=-l, -/? = 7 = 5, e"^" = . ^'^^^^ = ' ^^'^^ 

with (To 7^ 0. 

The same procedure can be followed to derive the singular solutions in the 
limit t — s> 0. The solution Aq is extended with the dilaton field going to a constant, 

$ = 00 + 01^^ 

a = l ,/9 = 7 = 2 

(^0): Sejao) 3e^° (4-10) 

52^'((Jo)2 + 3e2<^o ' ((jo)2 + 3e2<^o 



15 

27TU}f5 ' 

(4.5) 
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Furthermore, for the same form of oj and 5 < 0, two new solutions can be obtained 
with the moduhis behaving logarithmicaUy cr = ctq + /31nt and the dilaton being 
either of the form $ = 0o + 7 In ^ 

(Bo) : a = 1 , = -7 , = -ge--s-(/3)-o ^ (4.11) 

or of the form ^ = (jyit^ 

(Co): a = l, |/?| =2,7 = -2,01 = -v^exp(^^^^^^) . (4.12) 

Bo is a weak couphng solution, while Cq is a strong coupling one. 

Finally, for the same logarithmic behavior of lo and 5 > there is one more 
singular solution 

^ = 00 + 7lnt , 

(Do) : -(3 (4.13) 



cr = (To exp ^ 



with 



7 = 2, 3q;2 - 3^2 + 5a - 1 = , (3 = 5i"{Q)a^{a - 1) > 

2(3a - 1) (4.14) 



exp(0o) = 



3(1 — a)a^ 



which leads to a ~ 0.223, (3 ~ 0.0195, 0o ~ 3.24. This is a strong coupling 
asymptotic solution where the modulus field approaches the self-dual point in a 
non-analytic way. It provides an example of realizing the mechanism described in 
Section 3 according to which the Gauss-Bonnet term in (2.14) acts as a potential 
for the modulus with a minimum at the self-dual point. 

The integration of the non- linear system (2.12), (2.14), (2.15) can be performed 
numerically following a similar procedure we used in the previous Section for the 
metric-modulus case. However, in the presence of the dilaton the phase space is 
enlarged considerably since two more initial values are required for $ and $ in 
addition to a, uj and Co. Instead of presenting a detailed numerical investigation 
which will not be very illuminating, we concentrate to the analysis of some physically 
interesting cases. In the case 5 > 0, we verified that the fixed modulus singular 
solution Do (4.14) is obtained for a large region of the parameter space. It turns 
out that the modulus field approaches at t = its self-dual point by dumping 
oscillations consistently with the effective potential interpretation (see fig. 6). 
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The case 5 < is more interesting since it admits non-singular solutions. They 
could be in principle derived by the method described in Section 3 which consists 
of starting the numerical integration at a point where the energy conditions are 
violated, for instance when Co = 0. It turns out that this method is not very useful 
because, in this case, the violation of energy conditions although necessary is not 
sufficient to avoid the initial singularity [6]. An alternative way would be to start 
the integration with an asymptotically fiat solution in the infinite past, which in 
the metric-modulus case was shown to lead always to non-singular solutions. This 
turns out to be the case even in the presence of the dilaton. Starting from C^, 
which extends Cqo with a negligible dilaton, one is smoothly driven to A'^ in the 
infinite future. This confirms that the main characteristics of such solutions are a 
consequence of the modulus dependent string loop correction to the Gauss-Bonnet 
term and they do not depend on the existence of the dilaton. In a typical solution of 
this kind the scale factor and the internal radius behave similarly to those obtained 
in the absence of the dilaton (figs. 1,2). The dilaton evolution is presented in fig. 7; 
it starts from a constant value in the remote past and grows logarithmically towards 
strong coupling in the future. During the inflationary period when the modulus 
passes through its self-dual point, the dilaton also jumps to its maximum value. This 
can be easily understood by inspection of its equation of motion (2.15), where the 
Gauss-Bonnet term plays the role of a runaway potential for uj ~ constant. In the 
limit t — > cxD, although a weak coupling solution of the form A'^ or asymptotically 
constant dilaton of the form A^ are not a priori excluded, strong coupling seems to 
be preferred at least in the restricted region of the parameter space we scanned. In 
any case, at late times the effective action (2.6) should probably be modified by the 
addition of a dilaton and modulus potential arising from supersymmetry breaking 
or other non perturbative effects, which would stabilize these fields. 



5. Conclusions 



In this work we have examined the cosmological implications of the moduli 
dependent loop corrections to the gravitational couplings of the superstring effective 
action in the case of orbifold compactifications. These corrections consist of the 
Gauss-Bonnet integrand multiplied by a universal non-trivial function of the moduli 
fields and a numerical coefficient 5 which depends on the massless spectrum of 
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every particular model. We first derived the equations of motion for the metric, 
the dilaton and the modulus corresponding to the common compactification radius, 
and we classified all asymptotic solutions. Among them, there is one where the 
internal radius is fixed at its self-dual point as the universe approaches the initial 
singularity. 

In the case of negative sign for the parameter 6, we have shown that the strong 
energy condition p + p > associated to the stress energy tensor of the modulus 
can be violated, leading to an inflationary period and providing the possibility of 
avoiding the initial singularity. In fact a numerical analysis of the system has veri- 
fied the existence of a class of non-singular solutions which interpolate between an 
asymptotically fiat and a slowly expanding universe with a period of rapid expan- 
sion, when the modulus field passes through its self-dual value. This is in contrast 
to the behavior of other non-singular solutions which have been proposed in the 
literature, where de Sitter space was always obtained at "early" times [8,9], but it 
shares the properties of a "pre-big-bang" phase proposed in ref. [2]. 

Our solutions depend crucially on the form of the modulus dependent string 
loop corrections, while the dilaton contribution is negligible and can be ignored. 
Furthermore all time derivatives can remain bounded in the perturbative regime, 
which is consistent with our approximation of neglecting higher derivative terms in 
the effective action. Finally, the whole parameter range has been explored in the 
absence of the dilaton by the study of the corresponding phase diagram and we 
found that the class of non-singular solutions extends over an infinite region of the 
phase space. 
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Figure captions 

Fig. 1. The scale factor (continuous line) and the Huble expansion rate oj 
(dashed line) for the non-singular solution (5 < 0). 

Fig. 2. The modulus field a (continuous line) and its derivative a (dashed line) 
for the non-singular solution. 

Fig. 3. p + p (dashed line) and p + 3p (continuous line) for the non-singular solu- 
tion. 

Fig. 4. The maximum of the expansion rate cbmax as a function of the modulus a 

for o" > (continuous line) and d < (dashed line). 

Fig. 5. The phase diagram of cj as a function of a; continuous lines correspond to 
non-singular solutions and dashed lines to singular ones. The boundary of the two 

regions is plotted by a bold line. 

Fig. 6. The scale factor (dashed line) and the modulus a (continuous line) in a 
singular solution {5 > 0) where the modulus approaches its self-dual point (u = 0) 
at early times. 

Fig. 7. The dilaton field $ (continuous line) and its derivative $ (dashed line) for 
the non-singular solution. 
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